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We study the 0^"'" glueball correlator constructed with SU(3) anisotropic quenched lattice 
, QCD at various temperature taking into account the possible existence of the thermal width 

■ in the ground-state peak. For this purpose, we adopt the Breit-Wigner ansatz, analysing the 

lattice data obtained with 5,500-9,900 gauge configurations at each T. The results indicate 
the significant thermal width broadening as r(Tc) ~ 300 MeV with a reduction in the peak 
center as AuoiTc) ~ 100 MeV in the vicinity of the critical temperature Tc. 



§1. Introduction 



> 

o 

o , 

^ ■ At finite temperature/density, QCD changes its vacuum structure sucli as the 

\ reduction of the string tension, the partial chiral restoration, etc, even below the 

C3 ' critical temperature T^ of the QCD phase transition. These changes are expected to 

affect the hadrons leading to the changes in the various hadronic properties, since 
hadrons are composite particles consisting of quarks and gluons. The hadronic pole- 
^ ' mass shifts are thus considered to serve as the important precritical phenomena of 

^ . the QCD phase transition near the critical temperature Tc, and have been exten- 

sively studied by using various QCD-motivated low-energy effective theories.C^''^''^''^ 
These studies suggest the pole-mass reductions of charmoniums, light qq mesons and 
^ . the glueball near the critical temperature. 

^ \ As for the lattice QCD, it has been more popular to study the spatial correla- 

tions (the screening mass) than the temporal correlations (the pole-mass) at finite 
temperature. This is because the temporal extension of the lattice is subject to l/T. 
Accordingly, number of the temporal lattice data decreases at finite temperature, 
which used to be the main obstacle for the lattice QCD to measure the hadronic 
pole-mass. Recently, by using the anisotropic lattice, which has a finer lattice spac- 
ing in the temporal direction than in the spatial direction, the accurate measurement 
of the pole-mass with the lattice QCD has become possible at finite temperature.'^''^ 
Quenched-level Monte Carlo calculations reveal the profound results that the pole- 
masses of the qq mesons are almost unchanged from their zero-temperature values 
in the confinement phase, while the pole-mass of the O^"*" glueball shows the 300 
MeV reduction near the critical temperature.® Although these thermal effects are 
weaker than anticipated by the effective model studies, these tendencies are consis- 
tent with the recent quenched-level lattice studies on the screening mass at finite 
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temperatureP'^ In these analysis, the bound-state peaks in the spectral function 
are assumed to be narrow enough. However, at finite temperature, each bound-state 
acquires the thermal width through the interaction with the heat bath. The ther- 
mal width is expected to grow up with temperature, which may lead to a possible 
collapse of the narrow-peak assumption in some cases. In this paper, we first dis- 
cuss what is the expected consequence of the thermal width broadening. Then, we 
propose the Breit-Wigner ansatz for the fit-function for the temporal correlator at 
finite temperature. We finally show the results of the Breit-Wigner analysis of the 
0"*"+ glueball correlator at finite temperature.^ 



§2. The Breit-Wigner Ansatz 



We begin by considering the temporal correlator G(r) = Z{P) ^Tr (e ^^0(r)(/)(O)) 
and its spectral representation as 

^ .)) 

^ ' y_»2» 2siiih(/3ij/2) '^^ ^ ' 

where H denotes the QCD Hamiltonian, Z{I3) = Tr(e^^^) the partition func- 
tion, (/'(t) is the zero-momentum projected glueball operatoJ^ represented in the 
imaginary-time Heisenberg picture as (/'(r) = e'^^ (p{0)e~'^^ . Here, pico) denotes the 
spectral function 

Pi^) = E Kl^^^e-^^™ X 27r (5(a; - AE^m) - 8(0; + AEmn)) , (2-2) 

m,n ^ ' 

where En denotes the energy of nth excited states, and AEmn = Em — En- Note 
that p{io) is odd in u) reflecting the bosonic nature of the glueball. By adopting the 
appropriate ansatz for we can extract various physical quantities such as the 

pole-mass and the width through the spectral representation Eq. 1)2 -IJ) . 

We first consider the case where the bound-state peak is narrow .'^'^^''^ In this 
case, by introducing the temperature-dependent pole-mass m{T), p{uj) can be pa- 
rameterized as 

p{uj) ~ 2ttA {5{uj - m{T)) - 5{lo + m{T))) , (2-3) 

where A represents the strength. The second delta-function is introduced to respect 
the odd- function nature of p{ijj). Since the corresponding G(r) reduces to a single hy- 
perbolic cosine, the pole-mass measurement at finite temperature can be performed 
in the same way as the standard mass measurement at zero temperature.^ 

We next consider the case where the thermal width is wide. In this case, the 
peak center cjq of p{oj) represents the observed "mass" of the thermal hadron. What 
follows the narrow-peak assumption in this case ? To consider this, we notice that 
G(r) can be thought of as a weighted average of hyperbolic cosines with the weight 
as 
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Here, 2 sinh(/?u;/2) in the denominator works as a biased factor, which enhances the 
smaher uj region while suppressing the larger lo region. Consequently, the pole-mass 
m{T), which is approximated with the peak position of W{uj), is smaller than the 
peak center loq of the spectral function p{uj), i.e., the observed hadron "mass".^ 
What is the appropriate functional form of the fit-function ? To find this, we consider 
the retarded Green function Gr(u;). At T = 0, bound-state poles of GR(a;) are 
located on the real cj-axis. At T > 0, bound-state poles are moving into the complex 
(j-plane with increasing temperature. Suppose that a bound-state pole is located at 
uj = ujq — ir as 

GrH = ^—^ + --- , (2-5) 

where A represents the residue at the pole, and "• • • " the non-singular terms around 
the pole. Since the spectral function is the imaginary part of the retarded Green 
function, the contribution of this complex pole is expressed in the form of Lorentzian 

as 



piu;) = -2Im [Gniuj)) . (2-6) 
~ 2ttA (^5riio - loq) - 6riuJ + wo)) , 

where 6e{x) = — Im ( ) = is a smeared delta-function with the width 

vr \x — ie J ir x'' + 

e > 0. The second term in Eq. (|2-5() is introduced to respect the odd-function nature 

of p{uj). In the limit T ->■ +0, Eq. ((T^ reduces to Eq. (|^ . We thus see that the 

appropriate fit-function, which takes into account the effect of the non-zero thermal 

width, is the following Breit-Wigner type as 

, , f°° doj cosh(uj((3/2 - t)) ^ , , , x\ 



oo 



where A, F and are understood as fit-parameters, corresponding to the residue, 
the thermal width and the peak center, respectively. Note that g(r) is a generaliza- 
tion of the ordinary single hyperbolic cosine fit-function. In order to use Eas. H2-3() 
and ()2-6() . it is essential to suppress the higher spectral contributions. This is usually 
achieved by appropriate choices of the fit-ranges and also by improving the glueball 
operator, for instance, with the smearing method.!^ 



§3. Numerical Result 



We use the SU(3) anisotropic lattice plaquette action as 

^^ = W^ ^ ReTV(l-P,,(.)) + ^7GEl^«Tr(l-P,4(s)), (3-1) 

^'^ s,i<j<3 i,j<3 

where P^v{s) G SU(3) denotes the plaquette operator in the ;U-z/-plane. The lattice 
parameter and the bare anisotropic parameter are fixed as /3iat = '^^d 9^ = 6.25 and 
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7g = 3.2552, respectively, so as to reproduce renormalized anisotropy ^ = as /at = 
^pjnj These parameter set reproduces o^^ = 2.341(16) GeV and a^^ = 9.365(66) 
GeV, where the scale unit is introduced from the on-axis data of the static inter- 
quark potential with the string tension ^/a = 440 MeV. Numerical calculations are 
performed on the lattice of the size 20^ xNt with various A^^. The critical temperature 
Tc on this lattice is estimated from the behavior of the Polyakov-loop susceptibility 
as Tc ~ 280 MeV. The pseudo-heat-bath algorithm is adopted for the update of the 
gauge configurations. In order to construct the temporal glueball correlators, we use 
5,500 to 9,900 gauge configurations. The statistical data are divided into bins of 
the size 100 to reduce the possible auto-correlations near the critical temperature. 
The smearing method is used to obtain the improved glueball operator, which is 
determined by examining its behavior at the lowest temperature, i.e., T = 130 MeV. 
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Fig. 1. The center loq and the thermal width F of the lowest 0"'"'' glueball peak plotted against 
temperature T. The vertical dotted lines indicate the critical temperature Tc ~ 280 MeV. 

The peak center and the thermal width of thermal glueball are obtained by the 
best-fit analysis of the suitably smeared glueball correlator G(t)/G(0) with Eq. l\2-7\t 
at various temperatures. In Fig^ the peak center ujq and the thermal width F 
are plotted against temperature. While narrow-peak ansatz leads to the pole-mass 
reduction of 300 MeV near Tc in RelP the Breit-Wi gner analysis indicates a small 
reduction in the peak center as Aloo{Tc) ~ 100 MeV. Instead, we observe a significant 
thermal width broadening as r{Tc) ~ 300 MeV. 

§4. Summary and Discussion 

We have studied the temporal correlator of the 0"^~^ glueball at finite temperature 
using SU(3) anisotropic lattice QCD at quenched level with 5,500 to 9,900 gauge 
configurations at each temperature. We have proposed the Breit-Wigner ansatz for 
the fit-function to take into account the effect of the non-zero thermal width at finite 
temperature. We have applied the Breit-Wigner analysis to the temporal glueball 
correlator at finite temperature, and have observed a significant broadening of the 
thermal width as r{Tc) ~ 300 MeV with a slight reduction of the peak center as 
Au)q{Tc) ~ 100 MeV in the vicinity of the critical temperature Tc- 

In our present analysis, the Breit-Wigner ansatz seems to be the most appropri- 
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ate one. Above Tc, however, the particle spectrum is expected to change drastically, 
and a proper assumption on the spectral function becomes less trivial. In this re- 
spect, it is interesting to apply the recently developed maximum entropy method 
(MEM),E3ni which requires no assumption on the shape of the spectral function. 

Comparing our results^'' ' with the related lattice QCD results,'^''^ the quenched 
lattice QCD indicates that no significant thermal effects on the pole-mass have been 
so far observed in the qq sector,'^'^^ whereas a significant thermal effects can be ob- 
served in the glueball sector. These behaviors are consistent with the thermal effects 
on the screening mass.'^'^ 

To discuss the experimental observability, it is important to estimate the effects 
of the dynamical quarks. One of the most significant contributions of the light 
dynamical quarks in unquenched QCD is the "string breaking" phenomenon, i.e., the 
light dynamical quark screens the colored flux and the long flux tubes are expected to 
be broken into pieces. The string breaking thus affects the long distance behavior of 
the inter-particle interactions, which could lead to a possible changes in the structures 
of hadrons in quenched QCD. In the case of the O"*"^ glueball, since its size is known 
to be rather compact as -R ~ 0.4 fni,8'' -^^'' such changes of the glueball structure is 
expected to be less signiflcant. As for the mixing effect with the 0"^^ qq mesons 
in unquenched QCD, a recent unquenched lattice QCD study shows that there is 
no quark mass dependence on the glueball mass, which suggests that the mixing is 
actually rather small. However, it is in principle desirable to estimate the size of this 
mixing explicitly, which is left for the future studies. 
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